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Lecture 5 



In the previous lecture, the quantum mechanical origin of the linear and nonlinear susceptibilities 
was discussed. In particular, a perturbation analysis of the density operator was performed, and 
the resulting system of equations was solved recursively for the nth order density operator Pn{t) 
in terms of /5„_i(t), where the zeroth order term (independent of the applied electric field of the 
light) is given by the Boltzmann distribution at thermal equilibrium. 

So far we have obtained the linear optical properties of the medium, in terms of the first order 
susceptibility tensor (of rank- two), and we will now proceed with the next order of interaction, 
giving the second order electric susceptibility tensor (of rank three). 

The second order polarization density 

For the second order interaction, the corresponding term in the perturbation series of the density 
operator in the interaction picture becomes^ 



m = ^J jHi{n),p'i{Ti)]dn 



1 



[^i'(ri),[^i'(T2),Po]]rfr2dTi 



(Zfi)2 

In order to simplify the expression for the second order susceptibility, we will in the following 
analysis make use of a generalization of the cyclic perturbation of the terms in the commutator 
inside the trace, as 

Tr{[Qa(Tl),[g/3(T2),/5o]]g4 = Tr{po[[Q^,Qa(Ti)],Q;3(T2)]}. (2) 

By inserting the expression for the second order term of the perturbation series for the density 
operator into the quantum mechanical trace of the second order electric polarization density of the 
medium, one obtains 



Pj^\v,t) = ^TT[p2{t)Q, 



V 

iTr 
V 



{Uoit) ^ J p [^i'(n), [^i'(t2), po]] dT2 dn Uoi-t)) Q, 



=P2(t) (interaction picture) 



=p2(t) (Schrodinger picture) 
= {ii^<j('''i) and Ea{Ti) are classical fields (omit space dependence r)} 

= Tr {Uo{t) y* P' [g„(Ti), [Q0{t2), po]]Ea{ri)E0{T2) dT2 dnUo{-t) Q^] 

= {Pull out Eai{Ti)Ea2{'''2) and the integrals outside the trace} 



^ It should be noticed that the form given in Eq. (1) not only appUes to an ensemble of molecules, 
of arbitrary composition, but also to any kind of level of approximation for the interaction, such as 
the inclusion of magnetic dipolar interactions or electric quadrupolar interactions as well. These 
interactions should (of course) be incorporated in the expression for the interaction Hamiltonian 
H[{t), here described in the interaction picture. 
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Tr 

OO ^ — OO 



{ Uo (t) [Qc (n) , [Qp (t2) , Po]]M-t) (^l)^/3 (^2) dT2 ■ 



v{ihY 

In analogy with the results as obtained for the first order (linear) optical properties, now express 
the term -Bai (ti)£^q2 ("^2) in the frequency domain, by using the Fourier identity 



/oo 
-00 



which hence gives the second order polarization density expressed in terms of the electric field in 
the frequency domain as 

(^'*)= Wl^ r r f r Tr{?7o(t)[Q«(ri),[Q^(r2),/5o]][>oH)OM|£^a('^i)^^/3('^2) 
y ymj J -00 J -00 J -00 J —00 ^ ^ 

X exp(— iwiTi) exp(— ia;2T2) dT2 dri dw2 dwi 

= {Use exp(— iwr) = exp(— exp[— ia;(T — t)\} 

= I -00 /-I /-oo /-I {^°^*^[^«(^^^' /5o]]f/o(-t) Q^}E„{uJ^)Ep{uJ2) 

X exp[— ia;i(ri —t) — iu2{T2 — t)] dT2 dri exp[— i (wi + 0)2) t] duj2 divi 

^ ' 

/OO poo 
/ X^f^li3i-'^a;(^i,(^2)Ea{(Ji)Ei3{(j2) exp(-ja;<^t) duj2 dwi, 
-00 J — 00 

where the second order (quadratic) electric susceptibility is defined as 

(2) / \ 

= WiiKf /^/^^{^oW[^«(^i)'[^M^2),/5o]]J7o(-t)Q^} 

X exp[— ia;i(Ti —t) — iu)2{T2 — t)] dT2 dri 
= {Make use of Eq. (2) and take t{ = n — t} 



Tr{po[[Q^, Qairi)], Q^i^)]} exp[-j(a;iT{ + War^)] dr^ dT[. 



—00 J —00 



This obtained expression for the second order electric susceptibility does not possess the property 
of intrinsic permutation symmetry. However, by using the same arguments as discussed in the 
analysis of the polarization response functions in lecture two, we can easily verify that this tensor 
can be cast into a symmetric and antisymmetric part as 



^ ^ V ' 

symmetric part 

+ \ [xS/3(-'^<T;wi,a;2) -xlf/3a(-^<^;'^2,wi)], 
^ V " 



antisymmetric part 

and since the antisymmetric part, again following the arguments for the second order polarization 
response function, does not contribute to the polarization density, it is customary (in the Butcher 
and Cotter convention as well as all other conventions in nonlinear optics) to cast the second order 
susceptibility into the form 

xjfa/3(-t^a;Wl,a;2) 



ei,V{i 



j J ' Tr{po[[Q^,Q^{T[)],Qp{T^)]}exp[-i{uiT[+U2T^)]dT^dTi 
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where S, commonly called the symmetrizing operator, denotes that the expression that follows is 
to be summed over the 2! = 2 possible pairwise permutations of {ol^uji) and (/?, a;2), hence ensuring 
that the second order susceptibility possesses the intrinsic permutation symmetry, 

(2) / ^ (2) / \ 

X^a/3(-^<^;'^l''^2) = X^/3a(-^<^;'^2,Wi). 

Higher order polarization densities 

The previously described principle of deriving the susceptibilities of first and second order are 
straightforward to extend to the nth order interaction. In this case, we will make use of the 
following generalization of Eq. (2), 

Tr{[g„,(Ti), [Qa,(T2), . . . , [Qo„(r„),po]] • • 

= Tr{po[. . . [[Q^, Oai(ri)], Qa,(T2)], . . . Oa„(r„)]}, 

which, when applied in the evaluation of the expectation value of the electric dipole operator of 
the ensemble, gives the nth order electric susceptibility as 

(") { \ 

= 7W^^S /" r •■• r""Tr{po[...[[Q;.Qo,(ri)],Q„,(T2)],...Q„„(T„)]} 

£0<' \ '■''') J —oo J —oo J -oo 

X exp[-i(a;iTi + a;2T2 + . - . + WnT„)] drn ■■■ dT2 dri, 

where now the symmetrizing operator S indicates that the expression following it should be summed 
over all the n! pairwise permutations of (ai, Wi), . . . , (a„, a;„). 

It should be emphasized the symmetrizing operator S always implies summation over all the 
n! pairwise permutations of (ai, wi), . . . , (a„, a;„), irrcgardlcss of whether the permutations are 
distinct or not. This is due to that eventually occuring degenerate permutations are taken care of 
in the degeneracy coefficient K{—oj„;u}i, . . . , w„) in Butcher and Cotters convention, as described 
in lecture three and in additional notes that has been handed out in lecture four. 
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